Vector Calculus 2 



There's more to the subject of vector calculus than the material in chapter nine. There are a couple 
of types of line integrals and there are some basic theorems that relate the integrals to the derivatives, 
sort of like the fundamental theorem of calculus that relates the integral to the anti-derivative in one 
dimension. 

13.1 Integrals 

Recall the definition of the Riemann integral from section 1.6. 

/ dxf{x)= lim J2fi^k)^Xk (13.1) 

This refers to a function of a single variable, integrated along that one dimension. 

The basic idea is that you divide a complicated thing into little pieces to get an approximate 
answer. Then you refine the pieces into still smaller ones to improve the answer and finally take the 
limit as the approximation becomes perfect. 

What is the length of a curve in the plane? Divide the curve into a lot of small pieces, then if the 
pieces are small enough you can use the Pythagorean Theorem to estimate the length of each piece. 



A4 = ViAx^r + {^Vk? 



The whole curve then has a length that you estimate to be the sum of all these intervals. Finally take 
the limit to get the exact answer. 

5^ A4 = E V(Ax,)2 + (Ay,)2 I M = I Vdx^ + dy^ (13.2) 
k 

How do you actually do this? That will depend on the way that you use to describe the curve itself. 
Start with the simplest method and assume that you have a parametric representation of the curve: 

X = f{t) and y = g{t) 
Then dx = f{t)dt and dy = g{t)dt, so 



d^ = ^lUmf + {g{t)dtf = ^Jf{tf+g{tfdt (13.3) 
and the integral for the length is 

j dl = j\t^f{tY + g{tf 

where a and h are the limits on the parameter t. Think of this as J di = J v dt, where v is the speed. 
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Do the simplest example first. What is the circumference of a circle? Use the parametrization 

x = Rcos(j), y = RsiiKp then di = ^/ {-R sin (f)y + {R cos (pY d(j) = R dcj) (13.4) 

The circumference is then J dl = f^^ Rd(j) = 2'kR. An ellipse is a bit more of a challenge; see 
problem 13.3. 

If the curve is expressed in polar coordinates you may find another formulation preferable, though 
in essence it is the same. The Pythagorean Theorem is still applicable, but you have to see what it says 
in these coordinates. 



A4 = V(Arfc)2 + (rfcA0fc)2 




If this picture doesn't seem to show much of a right triangle, remember there's a limit involved, as Avf^ 
and A0/- approach zero this becomes more of a triangle. The integral for the length of a curve is then 



j M = j ^/(k^Tr^d^ 



To actually do this integral you will pick a parameter to represent the curve, and that parameter may 
even be (p itself. For an example, examine one loop of a logarithmic spiral; r = TqC^'^. 



di = Vrff2^T2d02 = ^[dr/d(t)Y + r"^ d(j) 
The length of the arc from = to = 27r is 

j ^J [roke^^f + (ro e^^f d(f) = j^"" d(j)ro e^^^/W+l = ro^Jk^ + l^ [e^^"" - l] 

If /c — )■ you can easily check that this give the correct answer. In the other extreme, for large k, you 
can also check that it is a plausible result, but it's a little harder to see. 

Weighted Integrals 

The time for a particle to travel along a short segment of a path is dt = dijv where v is the speed. 
The total time along a path is of course the integral of dt. 



dt 



;i3.5) 



How much time does it take a particle to slide down a curve under the influence of gravity? If the speed 
is determined by gravity without friction, you can use conservation of energy to compute the speed. I'll 
use the coordinate y measured downward from the top point of the curve, then 




mv"^ / 2 — mgy = E , so 



{2E/m) + 2gy 



(13.6) 
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Suppose that this particle starts at rest from y = 0, then E = and v = yJ2gy. Does the total time 
to reach a specific point depend on which path you take to get there? Very much so. This will lead to 
a classic problem called the "brachistochrone." See section 16.3 for that. 
1 Take the straight-line path from (0,0) to (xo,yo)- The path Is x = y ■ Xo/yo. 



V dx'^ + (i?/2 = dy^Jl+xl/y^, 



y^dy^l+xl/yl 



1 + xl/y^ 



1 1 



so 



2 There are an infinite number of possible paths, and another choice of path 
can give a smaller or a larger time. Take another path for which it's easy to 
compute the total time. Drop straight down in order to pick up speed, then V 
turn a sharp corner and coast horizontally. Compute the time along this path 
and it is the sum of two pieces. 



1 V^o + Vo 
2 



(13.7) 




y° dy dx _ 1 

Jo V^gm 



Xo 



1 



[xo + yo/2] 



(13.8) 



Which one takes a shorter time? See problem 13.9. 

3 What if the path is a parabola, x = y'^-Xo/y^? It drops rapidly at first, picking up speed, but then 
takes a more direct route to the end. Use y as the coordinate, then 

dx = 2y- Xo/yo, and di = .J{^ipxl/y^)^dy 



T 



dx 

V 



yo 



(4y 



o/yo 



+ 1 



-dy 



This is not an integral that you're likely to have encountered yet. I'll refer you to a large table of 
integrals, where you can perhaps find it under the heading of elliptic integrals. 

In more advanced treatments of optics, the time it takes light to travel along a path is of central 
importance because it is related to the phase of the light wave along that path. In that context however, 
you usually see it written with an extra factor of the speed of light. 



cT 



ndi 



(13.9) 



This last form, written in terms of the index of refraction, is called the optical path. Compare problems 
2.35 and 2.39. 

13.2 Line Integrals 

Work done on a point mass in one dimension is an integral. If the mass moves in three dimensions and 
if the force happens to be a constant, then work is simply a dot product: 



W 



Xf 



Fx{x) dx respectively W = F • Ar 



X; 



The general case for work on a particle moving along a trajectory in space is a line integral. It combines 
these two equations into a single expression for the work along an arbitrary path for an arbitrary force. 
There is not then any restriction to the elementary case of constant force. 
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The basic idea is a combination of Eqs. (13.1) and (13.2). Divide the specified curve into a 
number of pieces, at the points {r/j.}. Between points k — 1 and k you had the estimate of the arc 
length as ^/ (Ax^)^ + (Ayj^y, but here you need the whole vector from r^-i t° oi'der to evaluate 
the work done as the mass moves from one point to the next. Let Ar}, = r}, — Vj^^i, then 



^ ^ r ^ 

lim J]F(ffc)-Affc= Fir)' 



dr 




(13.10) 



This is the definition of a line integral. 

How do you evaluate these integrals? To repeat what happened with Eq. (13.2), that will depend 
on the way that you use to describe the curve itself. Start with the simplest method and assume that 
you have a parametric representation of the curve: r(t), then df= rdt and the integral is 



j F{r)-dr = j F{r{t))-fdt 



This is now an ordinary integral with respect to t. In many specific examples, you may find an easier 
way to represent the curve, but this is something that you can always fall back on. 

In order to see exactly where this is used, start with F = ma, Take the dot product with dr'and 
manipulate the expression. 



m 



dv 

~di'' 



so 
or 



F-dr 
F-dr 



dv ,^ 
m-j- ■ dr = 
dt 

—d{y-v) 



dv dr ,, ,^ dr ^ ,^ 
m-r- ■ -j-dt = mdv ■ -r = ttlv ■ dv 
dt dt dt 



(13.11) 



The integral of this from an initial point of the motion to a final point is 



F-df 



—d[v-v) = —[v^ 



(13.12) 



This is a standard form of the work-energy theorem in mechanics. In most cases you have to specify 
the whole path, not just the endpoints, so this way of writing the theorem is somewhat misleading. Is 
it legitimate to manipulate v • dv as in Eq. (13.11)? Yes. Simply write it in rectangular components as 
Vxdvx +Vydvy + Vzdvz and you can integrate each term with no problem; then assemble the result as 

Example 

If F = Axyx + B{x^ + L'^)y, what is the work done going from point (0, 0) to {L, L) along the three 
different paths indicated.? 



F-dr 



[FxX + Fyy] -[xdx + ydy] 

L rL 



F-dr 



F-df 



'C3 



[Fxdx + Fydy] = dxO+ dy 321"^ = 2BL^ 
Jo Jo 

[ dxAx^+f dyB{y^ + L^) = ALy3 + 4BLy3 
Jo Jo 

dyB{0 + L^)+ dxAxL = BL^ + AL^/2 
Jo Jo 
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Gradient 

What is the line integral of a gradient? Recall from section 8.5 and Eq. (8.16) that df = grad f ■ dr. 
The integral of the gradient is then 



^ grad/ ■ dr = j df = f2 - fi 



(13.13) 



where the indices represent the initial and final points. When you integrate a gradient, you need the 
function only at its endpoints. The path doesn't matter. Well, almost. See problem 13.19 for a caution. 

13.3 Gauss's Theorem 

The original definition of the divergence of a vector field is Eq. (9.9), 

div V = lim — —j— = lim — <t> v ■ dA 
V^o V dt v^o V J 

Fix a closed surface and evaluate the surface integral of v over that surface. 

dA 



v-dA 



S 




Now divide this volume into a lot of little volumes AV/^ with individual bounding surfaces Sf^. The 
picture on the right shows just two adjoining pieces of whole volume, but there are many more. If you 
do the surface integrals of v- dA over each of these pieces and add all of them, the result is the original 
surface integral. 



k 



V-dA 



v-dA 



s 



(13.14) 



The reason for this is that each interior face of volume V/j. is matched with the face of an adjoining 
volume Vi;/. The latter face will have dA pointing in the opposite direction, fi;^, = —fik, so when you 
add all the interior surface integrals they cancel. All that's left is the surface on the outside and the 
sum over all those faces is the original surface integral. 

In the equation (13.14) multiply and divide every term in the sum by the volume AV^. 



E 



V-dA 



Sk 



V-dA 



s 



Now increase the number of subdivisions, finally taking the limit as all the AVj. approach zero. The 
quantity inside the brackets becomes the definition of the divergence of v and you then get 



Gauss's Theorem: 



div V dV 



V 



V-dA 



(13.15) 



This* is Gauss's theorem, the divergence theorem. 



* You will sometimes see the notation dV instead of S for the boundary surface surrounding the 
volume V . Also dA instead of C for the boundary curve surrounding the area A. It is probably a 
better and more consistent notation, but it isn't yet as common in physics books. 
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Example 

Verify Gauss's Theorem for the solid hemisphere, r < R, 
< 6* < 7r/2, < < 27r. Use the vector field 

F = far'^sm9 + 9/3r9'^cos^(f) + 4)'jrsm9cos^(l) (13.16) 

Doing the surface integral on the hemisphere n = f, and on the bottom flat disk you have n = 9. The 
surface integral is then assembled from two pieces, 



(fF-dA= j far'^sm9-fdA+ j 9 l3r 9"^ cos'^ (f) ■ 9 dA 

J >/hemisph ^disk 

/■7r/2 1-2% pR p2-K 

= / R^sm9d9 / d^aR^sm9+ / rdr rf0/3r(7r/2)2 cos^ 
Jo Jo Jo Jo 

= a'K^R^/2 + /3'K^R^/l2 (13.17) 
Now do the volume integral of the divergence, using Eq. (9.16). 

divF = ^^ar^ sin 6' H —^-^P'r9'^ sm6' cos^ 6 H —p^Srl^ sin6'cos^ 6 

or r sin 9 o9 r sm 9 0(p 

= 4ar sm9 + /3 cos^ (f)[29 + 9'^ cot 9] + 27 sin cos 

The 7 term in the volume integral is zero because the 2 sin cos = sin 20 factor averages to zero 
over all 0. 

rR rTv/2 /■27r 

/ c/rrM sm9d9 dcp [Aar sin 9 + /3 cos^ (l)[29 + 9"^ cot 9]] 
Jo Jo Jo 

TT R3 fTT/2 

= 4a~-27i-j+l3-^-7vJ d9 sin 9[29 + 9'^ cot 9] 



aTi'^R^/2 + /Sn^Ryu 



The last integration used parametric differentiation starting from j^^'^ d9 cosk9, with differentiation 
with respect to k. 

13.4 Stokes' Theorem 

The expression for the curl in terms of integrals is Eq. (9.17), 

curl-u = lim ^ <L dA x v (13.18) 
V^o V J 

Use the same reasoning that leads from the definition of the divergence to Eqs. (13.14) and (13.15) 
(see problem 13.6), and this leads to the analog of Gauss's theorem, but with cross products. 



dAxv= / cnAvdV (13.19) 
S Jv 

This isn't yet in a form that is all that convenient, and a special case is both easier to interpret and 
more useful in applications. First apply it to a particular volume, one that is very thin and small. Take 
a tiny disk of height A/i, with top and bottom area AAi. Let fii be the unit normal vector out of the 
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top area. For small enough values of these dimensions, the right side of Eq. (13.18) is simply the value 
of the vector cmlv inside the volume times the volume AAiAh itself. 

(f dA X V = I cm\v dV = cm\v AAiAh f 

Js Jv 

Take the dot product of both sides with fii, and the parts of the surface integral from the top and 
the bottom faces disappear. That's just the statement that on the top and the bottom, dA is in the 
direction of ±ni, so the cross product makes dA x v perpendicular to -hi. 

I'm using the subscript i for the top surface and I'll use 2 for the surface around the edge. 
Otherwise it's too easy to get the notation mixed up. 

Now look at dA x v around the thin edge. The element of area has height Ah and length A£ 
along the arc. Call n2 the unit normal out of the edge. 



AA2 = AhAin2 

di 

The product ni ■ AA2 x v = ni ■ n2 x v AhAi = ni x n2-v AhAi, using the property of the triple 
scalar product. The product ni x n2 is in the direction along the arc of the edge, so 

nixn2Ai = Ai (13.20) 

Put all these pieces together and you have 



ni • d) dA X V = d) v-di Ah = fii ■ curl v AAiAh 

^ S c 

Divide by AAiAh and take the limit as AAi — )• 0. Recall that all the manipulations above work only 
under the assumption that you take this limit. 

fii ■ curl-u = lim -J—r(l)V-di (13.21) 
AA-^oAAJc 

You will sometimes see this equation (13.21) taken as the definition of the curl, and it does have an 
intuitive appeal. The one drawback to doing this is that it isn't at all obvious that the thing on the 
right-hand side is the dot product of Hi with anything. It is, because I deduced that fact from the 
vectors in Eq. (13.19), but if you use Eq. (13.21) as your starting point you have some proving to do. 

This form is easier to interpret than was the starting point with a volume integral. The line 
integral o^v-dd is called the circulation oft/ around the loop. Divide this by the area of the loop and 
take the limit as the area goes to zero and you then have the "circulation density" of the vector field. 
The component of the curl along some direction is then the circulation density around that direction. 
Notice that the equation (13.20) dictates the right-hand rule that the direction of integration around 
the loop is related to the direction of the normal fii. 

Stokes' theorem follows in a few lines from Eq. (13.21). Pick a surface A with a boundary C 
(or dA in the other notation). The surface doesn't have to be flat, but you have to be able to tell one 
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side from the other.* From here on I'll imitate the procedure of Eq. (13.14). Divide the surface into a 
lot of little pieces AAj^, and do the line integral v ■ d£ around each piece. Add all these pieces and 
the result is the whole line integral around the outside curve. 




(13.22) 



As before, on each interior boundary between area AAj^ and the adjoining AAj^/, the parts of the line 
integrals on the common boundary cancel because the directions of integration are opposite to each 
other. All that's left is the curve on the outside of the whole loop, and the sum over those intervals is 
the original line integral. 

Multiply and divide each term in the sum (13.22) by AA^. and you have 



E 



1 



AA 



k Jc. 



AA, 



v-di 



c 



(13.23) 



Now increase the number of subdivisions of the surface, finally taking the limit as all the AAi^ — )• 0, 
and the quantity inside the brackets becomes the normal component of the curl of v by Eq. (13.21). 
The limit of the sum is the definition of an integral, so 



Stokes' Theorem: 



curl V ■ dA 



A 



v-di 



c 



(13.24) 



What happens if the vector field v is the gradient of a function, v = V/? By Eq. (13.13) the 
line integral in (13.24) depends on just the endpoints of the path, but in this integral the initial and 
final points are the same. That makes the integral zero: fi — fi- That implies that the surface integral 
on the left is zero no matter what the surface spanning the contour is, and that can happen only if the 
thing being integrated is itself zero, curlgrad/ = 0. That's one of the common vector identities in 
problem 9.36. Of course this statement requires the usual assumption that there are no singularities of 
V within the area. 



Example 

Verify Stokes' 

< <i < 27r. 



theorem for that part of a spherical surface r 
Use for this example the vector field 



F = fAr"^ sin 9 + OBrO"^ cos (j) + 4>Crsin9 cos^ 

— * 

To compute the curl of F, use Eq. (9.33), getting 



R,0<9<9o 



(13.25) 




V X F = f 



1 



d 



r sin ^ V d6 



sin 6 Cr sin 6 cos^ i 



d 



cos ( 



+ 



. 1 

r sin^^ 



{Cr cos^ (j)2sm9 cos 9 + Br 9"^ sin 0) + 



(13.26) 



That means no Klein bottles or Mobius strips. 
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I need only the f component of the curl because the surface integral uses only the normal (f) 
component. The surface integral of this has the area element dA = sm.9 dO d(j). 



cuil F-dA= R^sinOde d<j)- ^(CRcos"^ (p2sm9 cos9 + BR9^ sine 

Jo Jo Rsm9^ 



/•Co /■27r 

= i?M d9 d(j)2Ccos'^(j)sm9cos9 
Jo Jo 

= R^2Ct: sin^ 9q/2 = CR\ sin^ 9o 
The other side of Stokes' theorem is the line integral around the circle at angle ^o- 



j)F-di= J rsm9od(j)Crsm9cos^ 



2tt 

■"'■CR^ sin^ ^ocos2, 







= CR' sm'9o7r (13.27) 

and the two sides of the theorem agree. Check! Did I get the overall signs right? The direction of 
integration around the loop matters. A further check: If = tt, the length of the loop is zero and 
both integrals give zero as they should. 

Conservative Fields 

An immediate corollary of Stokes' theorem is that if the curl of a vector field is zero throughout a region 
then line integrals are independent of path in that region. To state it a bit more precisely, in a volume 
for which any closed path can be shrunk to a point without leaving the region, if the curl of v equals 
zero, then F ■ df^ depends on the endpoints of the path, and not on how you get there. 
To see why this follows, take two integrals from point a to point b. 



V ■ df and v-df 
1 J2 



The difference of these two integrals is 

J V ■ df — J V ■ df= ^ V ■ df 

This equations happens because the minus sign is the same thing that you get by integrating in the 
reverse direction. For a field with V x v = 0, Stokes' theorem says that this closed path integral is 
zero, and the statement is proved. 

What was that fussy-sounding statement "for which any closed path can be shrunk to a point 
without leaving the region?" Consider the vector field in three dimensions, written in rectangular and 
cylindrical coordinates, 

V = A{xy-yx)/{x^ + y^) = A4>/r (13.28) 

You can verify (in either coordinate system) that its curl is zero — except for the ^-axis, where it is 
singular. A closed loop line integral that doesn't encircle the 2;-axis will be zero, but if it does go around 
the axis then it is not. The vector's direction 9 always points counterclockwise around the axis. See 
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problem 13.17. If you have a loop that encloses the singular line, then you can't shrink the loop without 
its getting hung up on the axis. 

The converse of this theorem is also true. If every closed-path line integral of v is zero, and if 
the derivatives o\ v are continuous, then its curl is zero. Stokes' theorem tells you that every surface 
integral of V x is zero, so you can pick a point and a small /S.A at this point. For small enough 
area whatever the curl is, it won't change much. The integral over this small area is then V x iT- /S.A, 
and by assumption this is zero. It's zero for all values of the area vector. The only vector whose dot 
product with all vectors is zero is itself the zero vector. 

Potentials 

The relation between the vanishing curl and the fact that the line integral is independent of path leads 
to the existence of potential functions. 

If curli^ = in a simply-connected domain (that's one for which any closed loop can be shrunk 
to a point), then I can write as a gradient, — grad/. The minus sign is conventional. I've already 
constructed the answer (almost), and to complete the calculation note that line integrals are independent 
of path in such a domain, and that means that the integral 



r 



F-dr (13.29) 



ro 



is a function of the two endpoints alone. Fix and treat this as a function of the upper limit r. Cal 
it —f{f). The defining equation for the gradient is Eq. (8.16), 

df = grad / ■ dr 

How does the integral (13.29) change when you change fa bit? 



/ F-df- / F-df= / F-df=F-df 

Jfo Jfo Jr 



This is ~df because I called this integral —f{f). Compare the last two equations and because df \s 
arbitrary you immediately get 

F = -grad/ (13.30) 

I used this equation in section 9.9, stating that the existence of the gravitational potential energy 
followed from the fact that V x ^ = 0. 

Vector Potentials 

This is not strictly under the subject of conservative fields, but it's a convenient place to discuss it 
anyway. When a vector field has zero curl then it's a gradient. When a vector field has zero divergence 
then it's a curl. In both cases the converse is simple, and it's what you see first: V x V/ = and 
V- V X A = (problem 9.36). In Eqs. (13.29) and (13.30) I was able to construct the function / 
because V x = 0. It is also possible, if V ■ = 0, to construct the function A such that F = V x A. 

In both cases, there are extra conditions needed for the statements to be completely true. To 
conclude that a conservative field (V x = 0) is a gradient requires that the domain be simply- 
connected, allowing the line integral to be completely independent of path. To conclude that a field 
satisfying V ■ -F = can be written as = V x A requires something similar: that all closed surfaces 
can be shrunk to a point. This statement is not so easy to prove, and the explicit construction of A 
from F is not very enlightening. 
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You can easily verify that A = S x r/2 is a vector potential for the uniform field B. Neither the 
scalar potential nor the vector potential are unique. You can always add a constant to a scalar potential 
because the gradient of a scalar is zero and it doesn't change the result. For the vector potential you 
can add the gradient of an arbitrary function because that doesn't change the curl. 

F = -V(/ + C) = -V/, and B = V x {A + V f) = V x A (13.31) 



13.5 Reynolds Transport Theorem 

When an integral has limits that are functions of time, how do you differentiate it? That's pretty easy 
for one-dimensional integrals, as in Eqs. (1.19) and (1.21). 

jj^^^^ dxg{x,t) = j^^^^ dx^^^+g{f2{t),t)^-g{Mt),t)^ (13.32) 

One of Maxwell's equations for electromagnetism is 

dB 

VxE = -^ (13.33) 
dt 

Integrate this equation over the surface S. 

J VxE-dA = j E-di = J -^-dA (13.34) 

This used Stokes' theorem, and I would like to be able to pull the time derivative out of the integral, 
but can 1? If the surface is itself time independent then the answer is yes, but what if it isn't? What if 
the surface integral has a surface that is moving? Can this happen? That's how generators works, and 
you wouldn't be reading this now without the power they provide. The copper wire loops are rotating 
at high speed, and it is this motion that provides the EMF. 

I'll work backwards and compute the time derivative of a surface integral, allowing the surface 
itself to move. To do this, I'll return to the definition of a derivative. The time variable appears in two 
places, so use the standard trick of adding and subtracting a term, just as in section 1.5. Call $ the 
flux integral, J B ■ dA. 



A$= / B{t + At)-dA- [ B{t)-dA 

JS{t+At) Js(t) 

= [ B{t + At)-dA- [ B{t)-dA (13.35) 

JS{t+At) JS{t+At) 

+ / B{t)-dA- j B{t) 

JS(t+At) JS(t) 



■dA 

IS{t+At) JS{t) 

B '\s a function of too, but I won't write it. The first two terms have the same surface, so they 
combine to give 

/ AB-dA 

Jsit+At) 

and when you divide by At and let it approach zero, you get 



Js[t) dt 
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Now for the next two terms, which require some manipulation. Add and subtract the surface that forms 
the edge between the boundaries C{t) and C{t + At). 









Isi 



B{t)-dA- [ B{t)-dA= (( B{t)-dA- [ B-dA (13.36) 

ISit+At) Js{t) J J edge 

The strip around the edge between the two surfaces make the surface integral closed, but I then have 
to subtract it as a separate term. 

You can convert the surface integral to a volume integral with Gauss's theorem, but it's still 
necessary to figure out how to write the volume element. [Yes, V-B = 0, but this result can be 
applied in other cases too, so don't use that fact here.] The surface is moving at velocity v, so an area 
element AA will in time At sweep out a volume AA-vAt. Note: v isn't necessarily a constant in 
space and these surfaces aren't necessarily flat. 

AV = AA-vAt^ f B{t)-dA= [ d^rV-B = [ V-BdA-vAt (13.37) 

J J Js(t) 

To do the surface integral around the edge, use the same method as in deriving Stokes' theorem, 
Eq. (13.20). 



AA = A£x vAt ^-^^^SU^^Z 





f B-dA= I B-dIxvAt= j vxB-diAt (13.38) 

J edge J C J C 

Put Eqs. (13.37) and (13.38) into Eq. (13.36) and then into Eq. (13.35), divide by At and let At 0. 

^ / B-dA= [ ^■dA+ [ V-Bv-dA- [ vxB-di (13.39) 
dt Js{t) Js{t) <Ji Js{t) Jc(t) 

This transport theorem is the analog of Eq. (13.32) for a surface integral. 

In order to check this equation, and to see what the terms do, try some example vector functions 
that isolate the terms, so that only one of the terms on the right side of (13.39) is non-zero at a time. 

1: B = Bozt, with a surface z = 0, + < R"^ 

For a constant Bq, and v = 0, only the first term is present. The equation is BqTcR^ = BqT^E? . 
Now take a static field 

2 : B = Czz, with a moving surface z = vt, + < E? 

The first and third terms on the right vanish, and V • B = C . The other terms are 



^Czz ■ -nR^ z 
dt 



z=vt 



CvttR^ = J {C)vz ■ dA = CvttR^ 
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Now take a uniform static field 



3 : 



B = BqZ with a radially expanding surface z = 0, x^ + ip'KR^, 



R = vt 



The first and second terms on the right are now zero, and 



di 



{vf X Bo 2) -Ode 



{-vBoe)-edi = +vBo27rR 



R=vt 



2Bo7ivH 



Draw some pictures of these three cases to see if the pictures agree with the algebra. 
Faraday's Law 

If you now apply the transport theorem (13.39) to Maxwell's equation (13.34), and use the fact that 



This is Faraday's law, saying that the force per charge integrated around a closed loop (called the EMF) 
is the negative time derivative of the magnetic flux through the loop. 

Occasionally you will find an introductory physics text that writes Faraday's law without the 
V X B term. That's o.k. as long as the integrals involve only stationary curves and surfaces, but some 
will try to apply it to generators, with moving conductors. This results in amazing contortions to try 
to explain the results. For another of Maxwell's equations, see problem 13.30. 

The electromagnetic force on a charge \s F = q(^E + v x B^ This means that if a charge inside 
a conductor is free to move, the force on it comes from both the electric and the magnetic fields in this 
equation. (The Lorentz force law.) The integral of this force ■ di is the work done on a charge along 
some specified path. If this integral is independent of path: V x E = and v = 0, then this work 
divided by the charge is the potential difference, the voltage, between the initial and final points. In the 
more general case, where one or the other of these requirements is false, then it's given the somewhat 
antiquated name EMF, for "electromotive force." (It is often called "voltage" anyway, though if you're 
being fussy that's not really correct.) 

13.6 Fields as Vector Spaces 

It's sometimes useful to look back at the general idea of a vector space and to rephrase some common 
ideas in that language. Vector fields, such as E{x,y,z) can be added and multiplied by scalars. They 
form vector spaces, infinite dimensional of course. They even have a natural scalar product 



Here I'm assuming that the scalars are real numbers, though you can change that if you like. For this 
to make sense, you have to assume that the fields are square integrable, but for the case of electric or 
magnetic fields that just means that the total energy in the field is finite. Because these are supposed 
to satisfy some differential equations (Maxwell's), the derivative must also be square integrable, and I'll 
require that they go to zero at infinity faster than l/r'^ or so. 

The curl is an operator on this space, taking a vector field into another vector field. Recall the 
definitions of symmetric and hermitian operators from section 7.14. The curl satisfies the identity 



V ■ 5 = you get 




(13.40) 




(13.41) 



{Ei,VxE2) = {VxEi,E2) 



(13.42) 
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For a proof, just write it out and then find the vector identity that will allow you to integrate by parts. 

V ■ {Ax B) = B -V X A- A-V X B (13.43) 

Equation (13.42) is 

J d^r Ei{f)-V X E2{r) = j d^r [V x Ei{r)) ■ E2{r) - j d^rV-{Ei x E2) 

The last integral becomes a surface integral by Gauss's theorem, § dA - {Ei x E2), and you can now 
let the volume (and so the surface) go to infinity. The fields go to zero sufficiently fast, so this is zero 
and the result is proved: Curl is a symmetric operator. Its eigenvalues are real and its eigenvectors are 
orthogonal. This is not a result you will use often, but the next one is important. 

Helmholtz Decomposition 

There are subspaces in this vector space of fields: (1) The set of all fields that are gradients. (2) The 
set of all fields that are curls. These subspaces are orthogonal to each other; every vector in the first 
is orthogonal to every vector in the second. To prove this, just use the same vector identity (13.43) 
and let A = V/. I will first present a restricted version of this theorem because it's simpler. Assume 
that the domain is all space and that the fields and their derivatives all go to zero infinitely far away. 
A generalization to finite boundaries will be mentioned at the end. 

V/ ■ V X = ■ V X V/ - V ■ (V/ X 5) 
Calculate the scalar product of one vector field with the other. 

(V/, V xB) = J d^rVf -V X B = j d\ [B - V xV f ~V ■ [V f x B)] 

= 0- ^ {Vf X B) ■dA = (13.44) 

As usual, the boundary condition that the fields and their derivatives go to zero rapidly at infinity kills 
the surface term. This proves the result, that the two subspaces are mutually orthogonal. 

Do these two cases exhaust all possible vector fields? In this restricted case with no boundaries 
short of infinity, the answer is yes. The general case later will add other possibilities. Here you have 
two orthogonal subspaces, and to show that these two fill out the whole vector space, I will ask the 
question: what are all the vector fields orthogonal to both of them? I will show first that whatever they 
are will satisfy Laplace's equation, and then the fact that the fields go to zero at infinity will be enough 
to show that this third case is identically zero. This statement is the Helmholtz theorem: Such vector 
fields can be written as the sum of two orthogonal fields: a gradient, and a curl. 

To prove it, my plan of attack is to show that if a field F is orthogonal to all gradients and to all 
curls, then V'^F is orthogonal to all square-integrable vector fields. The only vector that is orthogonal 
to everything is the zero vector, so F satisfies Laplace's equation. The assumption now is that for 
general / and v, 

Jd^rF-Vf = and Jd^rF-Vxv = 



I want to show that for a general vector field u, 

d^r u ■ V^F = 



gradie^ats 
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The method is essentially two partial integrals, moving two derivatives from F over to u. Start with 
the d'^/dz^ term in the Laplacian and hold off on the dx and dy integrals. Remember that all these 
functions go to zero at infinity. Pick the i-component of u and the j-component of F. 



Q2 

dz U{ Q^2 ~ ^i^zFj 



dz {dzUi) {dzFj 



0^{d,u,)Fj 



/POO 
dz[d'lui)Fj = / dz[d'lui)Fj 
J —oo 



Now reinsert the dx and dy integrals. Repeat this for the other two terms in the Laplacian, d'^Fj and 



dyFj. The result is 



/ 



d^r u ■ V^F 



d^r (V'^u)-F 



(13.45) 



If this looks familiar it is just the three dimensional version of the manipulations that led to Eq. (5.15). 
Now use the identity 

V X (V X m) = V(V-m) - V^m 



in the right side of (13.45) to get 



(fru-V'^F 



d^r 



(V(V ■ n )) ■ F - (V X (V X n )) ■ F 



(13.46) 



The first term on the right is the scalar product of the vector field F with a gradient. The second term is 
the scalar product with a curl. Both are zero by the hypotheses of the theorem, thereby demonstrating 
that the Laplacian of F is orthogonal to everything, and so V^F = 0. 

When you do this in all of space, with the boundary conditions that the fields all go to zero at 
infinity, the only solutions to Laplace's equation are identically zero. In other words, the two vector 
spaces (the gradients and the curls) exhaust all the possibilities. How to prove this? Just pick a 
component, say Fx, treat it as simply a scalar function — call it / — and apply a vector identity, 
problem 9.36. 

V-(0i') = (V0)-i+0(V-i) 

Let0 = / and A = Vf, then v .(fy f) = y f .y f + fy^ f 
Integrate this over all space and apply Gauss's theorem. {I.e. integrate by parts.) 

I d\ V -(/V/) = j>fyf.dA = j d\ [V/ ■ V/ + /VV] (13.47) 

If / and its derivative go to zero fast enough at infinity (a modest requirement), the surface term, 
§ dA, goes to zero. The Laplacian term, V^/ = 0, and all that's left is 



j d^rVf-Vf = 



This is the integral of a quantity that can never be negative. The only way that the integral can be 
zero is that the integrand is zero. If V/ = 0, then / is a constant, and if it must also go to zero far 
away then that constant is zero. 
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This combination of results, the Helmholtz theorem, describes a field as the sum of a gradient 
and a curl, but is there a way to find these two components explicitly? Yes. 

F = V/ + VxS, so V-F = VV, and V x F = V x V x S = V(V ■ S ) - V^S 
Solutions of these equations are 

/(f) = Zl / d'r' ^:^^f? and B{f) = — [ dh' Z^Zg!! (13.48) 
■'^ ^ An J \t^-t^'\ ^ ^ An J \f-f'\ ^ ' 

Generalization 

In all this derivation, I assumed that the domain is all of three-dimensional space, and this made the 
calculations easier. A more general result lets you specify boundary conditions on some finite boundary 
and then a general vector field is the sum of as many as five classes of vector functions. This is the 
Helmholtz-Hodge decomposition theorem, and it has applications in the more complicated aspects of 
fluid flow (as if there are any simple ones), even in setting up techniques of numerical analysis for such 
problems. The details are involved, and I will simply refer you to a good review article* on the subject. 



Exercises 

— * — * — * 

1 For a circle, from the definition of the integral, what is § Ml What is § Ml What \s § MxC where 

— * 

C is a constant vector? 

2 What is the work you must do in lifting a mass m in the Earth's gravitational field from a radius Ri 
to a radius R2- These are measured from the center of the Earth and the motion is purely radial. 

3 Same as the preceding exercise but the motion is 1. due north a distance Ri9q then 2. radially out 
to R2 then 3. due south a distance R2O0. 

4 Verify Stokes' Theorem by separately calculating the left and the right sides of the theorem for the 
case of the vector field 

— * 

F{x,y) = xAy + yBx 
around the rectangle {a < x < b), {c < y < d). 

5 Verify Stokes' Theorem by separately calculating the left and the right sides of the theorem for the 
case of the vector field 

F{x, y) = xAy — yBx 
around the rectangle (a < x <b), {c < y < d). 

6 Verify Stokes' Theorem for the semi-cylinder Q < z < h, < (f) < n, r = R. The vector field is 

F{r, 0, z) = fAr'^ sin + 4>Br(j)'^z + zCrz'^ 

7 Verify Gauss's Theorem using the whole cylinder Q < z < h, r = R and the vector field F{r, 0, z) = 

fAr"^ sin + (pBrz sin^ + zCrz"^ . 

8 What would happen if you used the volume of the preceding exercise and the field of the exercise 
before that one to check Gauss's law? 

* Cantarella, DeTurck, and Gluck: The American Mathematical Monthly, May 2002. The paper 
is an unusual mix of abstract topological methods and very concrete examples. It thereby gives you a 
fighting chance at the subject. 
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Problems 



13.1 In the equation (13.4) what happens if you start with a different parametrization for x and y, 
perhaps x = i?cos(0'/2) and y = i?sin(0'/2) for < (f)' < An. Do you get the same answer? 

13.2 What is the length of the arc of the parabola y = {o? — (—a < x < a)l 

But First draw a sketch and make a rough estimate of what the result ought to be. Then do the 
calculation and compare the answers. What limiting cases allow you to check your result? 
Ans: (6/2) [ sinh"^ c + cVl+~(?\ where c = 2a/h 

13.3 (a) You can describe an ellipse as x = acos0, y = bsm(f). (Prove this.) 

(b) Warm up by computing the area of the ellipse. 

(c) What is the circumference of this ellipse? You will find a (probably) unfamiliar integral here, so to 

put this integral into a standard form, note that it is 4 /J^^^. Then use cos^ = 1 — sin^0. Finally, 
look up chapter 17, Elliptic Integrals, of Abramowitz and Stegun. You will find the reference to this 
at the end of section 1.4. Notice in this integral that when you integrate, it will not matter whether 
you have a sin^ or a cos^. Ans: 4aE{m) 

13.4 For another derivation of the work-energy theorem, one that doesn't use the manipulations of 
calculus as in Eq. (13.11), go back to basics. 

(a) For a constant force, start from F = ma and derive by elementary manipulations that 



All that you need to do is to note that the acceleration is a constant so you can get v and r as functions 
of time. Then eliminate t 

(b) Along a specified curve Divide the curve at points 



In each of these intervals apply the preceding equation. This makes sense in that if the interval is small 
the force won't change much in the interval. 

(c) Add all these N equations and watch the kinetic energy terms telescope and (mostly) cancel. This 
limit as all the Af/, is Eq. (13.12). 

13.5 The manipulation in the final step of Eq. (13.12) seems almost too obvious. Is it? Well yes, but 
write out the definition of this integral as the limit of a sum to verify that it really is easy. 

13.6 Mimic the derivation of Gauss's theorem, Eq. (13.15), and derive the identities 



13.7 The force by a magnetic field on a small piece of wire, length di, and carrying a current / is 
dF = I di X B. The total force on a wire carrying this current in a complete circuit is the integral of 
this. Let B = X Ay — y Ax. The wire consists of the line segments around the rectangle < a; < a, 
< y < b. The direction of the current is in the +y direction on the x = line. What is the total 
force on the loop? Ans: 



F ■ Af = y [vf - vf] 



ri = ro, ri, r2 




and 
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13.8 Verify Stokes' theorem for the field F = Axyx + B{1 + x'^y'^)y and for the rectangular loop 

a < X <b, c < y < d. 

13.9 Which of the two times in Eqs. (13.7) and (13.8) is shorter. (Compare their squares; it's easier.) 

13.10 Write the equations (9.36) in an integral form. 

13.11 Start with Stokes' theorem and shrink the boundary curve to a point. That 
doesn't mean there's no surface left; it's not flat, remember. The surface is pinched off 
like a balloon. It is now a closed surface, and what is the value of this integral? Now 
apply Gauss's theorem to it and what do you get? Ans: See Eq. (9.34) 

13.12 Use the same surface as in the example, Eq. (13.25), and verify Stokes' theorem for the vector 
field 

F = f Ar~^ cos^ 9 sin (p + 6 Br^ sin 9 cos^ (p + 4>Cr~'^ cos^ 9 sin^ (f) 

13.13 Use the same surface as in the example, Eq. (13.25), and examine Stokes' theorem for the vector 
field 

F = ff{r, 9, 0) + 9g{r, 9, 0) + 4>h{r, 9, 0) 

(a) Show from the line integral part that the answer can depend only on the function h, not / or g. 

(b) Now examine the surface integral over this cap and show the same thing. 

13.14 For the vector field in the x-y plane; F = [xy — yx^ /2, use Stokes' theorem . — 

to compute the line integral of i^-rfr around an arbitrary closed curve. What is the \_y^ 
significance of the sign of the result? When you considered an "arbitrary" loop, did 
you consider the possibilities presented by these curves? 

13.15 What is the (closed) surface integral of = r/3 over an arbitrary closed surface? Ans; V . 

13.16 What is the (closed) surface integral of = r/3 over an arbitrary closed surface? This time 

however, the surface integral uses the cross product; § dA x F. If in doubt, try drawing the picture 
for a special case first. 

13.17 For the vector field Eq. (13.28) explicitly show that § v ■ dr is zero for a curve such 
as that in the picture and that it is not zero for a circle going around the singularity. 





13.18 Refer to Eq. (13.27) and check it for small 9q. Notice the combination tt{R9o 



|2 



13.19 For the vector field, Eq. (13.28), use Eq. (13.29) to try to construct a potential function. Because 
within a certain domain the integral is independent of path, you can pick the most convenient possible 
path, the one that makes the integration easiest. What goes wrong? 

13.20 Refer to problem 9.33 and construct the solutions by integration, using the methods of this 
chapter. 

13.21 (a) Evaluate § F ■ df ^or F = xAxy + yBx around the circle of radius R centered at the 
origin. 

(b) Now do it again, using Stokes' theorem this time. 

13.22 Same as the preceding problem, but § df x F instead. 
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13.23 Use the same field as the preceding two problems and evaluate the surface integral of F ■ dA 
over the hemispherical surface x'^ + y'^ + z"^ = E?, z > 0. 

13.24 The same field and surface as the preceding problem, but now the surface integral dA x F . 
Ans: z2nBr^/3 

13.25 (a) Prove the identity V ■ {A x B) = B -V x A - A -V x B . (index mechanics?) 

(b) Next apply Gauss's theorem to V ■ (A x i?) and take the special case that B is an arbitrary constant 
to derive Eq. (13.19). 

13.26 (a) Prove the identity V -(/F) = fV -F + F-Vf. 

(b) Apply Gauss's theorem to V-(/-F) for an arbitrary constant F to derive a result found in another 
problem. 

(c) Explain why the word "arbitrary" is necessary here. 

13.27 The vector potential is not unique, as you can add an arbitrary gradient to it without affecting 
its curl. Suppose that B = V x A with 

A = xaxyz + yf3x'^z + Z'jxyz'^ 

Find a function f{x,y,z) such that A' = A + Vf has the ^-component identically zero. Do you get 
the same B by taking the curl of A and of A'7 

13.28 Take the vector field 

B = axyx + f3xyy + ^{xz + yz)z 

Write out the equation B = 'VxA\n rectangular components and figure out what functions Ax{x, y, z), 
Ay{x,y, z), and Az{x,y,z) will work. Note: From the preceding problem you see that you may if you 

wish pick any one of the components of A to be zero — that will cut down on the labor. Also, you 
should expect that this problem is impossible unless B has zero divergence. That fact should come out 
of your calculations, so don't put it in yet. Determine the conditions on a, (3, and 7 that make this 
problem solvable, and show that this is equivalent to V ■ i? = 0. 

13.29 A magnetic monopole, if it exists, will have a magnetic field ^^q^f / Aur"^ . The divergence of this 
magnetic field is zero except at the origin, but that means that not every closed surface can be shrunk 
to a point without running into the singularity. The necessary condition for having a vector potential is 
not satisfied. Try to construct such a potential anyway. Assume a solution in spherical coordinates of 
the form A = <t> f{r)g{6) and figure out what / and g will have this B for a curl. Sketch the resulting 
A. You will run into a singularity (or two, depending). Ans: A = 0yUoQ'm(l — cos ^^)/(47rr^ sin0) (not 
unique) 

13.30 Apply the Reynolds transport theorem to the other of Maxwell's equations. 

V X B = /ioj +/ioeo-^ 

Don't simply leave the result in the first form that you find. Manipulate it into what seems to be the 
best form. Use yUo^o = l/c^- 

Ans: § [B - V X E/c^) ■di = fiQ J - pv) ■ dA + fioeo{d / dt) J E-dA 
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13.31 Derive the analog of the Reynolds transport theorem, Eq. (13.39), for a line integral around a 
closed loop. 

(a) ^/ F{f,t)-di= (f ^-di-f vx{VxF)-di 
dt Jc(t) Jc(t) Jc(t) 

and for the surface integral of a scalar. You will need problem 13.6. 

(b) ^/ (p{r,t)dA= [ ^dA+ [ {V (j)) dA ■ v - l (pdixv 

Js{t) Js{t) Ot Js{t) Jc{t) 

Make up examples that test the validity of individual terms in the equations. I recommend cylindrical 
coordinates for your examples. 



13.32 Another transport theorem is more difficult to derive. 




I had to look up some vector identities, including one for V x (A x i?). A trick that I found helpful: 
At a certain point take the dot product of the whole equation with a fixed vector B and manipulate 
the resulting product, finally factoring out the arbitrary vector B ■ at the end. Make up examples that 
test the validity of individual terms in the equations. Again, I recommend cylindrical coordinates for 
your examples. 

13.33 Apply Eq. (13.39) to the velocity field itself. That is, let B = v. Suppose further that the fluid 
is incompressible with V-{/ = and that the flow is stationary (no time dependence). Explain the 
results. 

13.34 Assume that the Earth's atmosphere obeys the density equation p = pQe~^/^ for a height z 
above the surface, (a) Through what amount of air does sunlight have to travel when coming from 
straight overhead? Take the measure of this to be f pdi (called the "air mass"), (b) Through what 
amount of air does sunlight have to travel when coming from just on the horizon at sunset? Neglect 
the fact that light will refract in the atmosphere and that the path in the second case won't really be a 
straight line. Take = 10 km and the radius of the Earth to be 6400 km. The integral you get for the 
second case is probably not familiar. You may evaluate it numerically for the numbers that I stated, 
or you may look it up in a big table of integrals such as Gradshteyn and Ryzhik, or you may use an 
approximation, <C -R. (I recommend the last.) What is the numerical value of the ratio of these two 
air mass integrals? This goes far in explaining why you can look at the setting sun. 

(c) If refraction in the atmosphere is included, the light will bend and pass through a still larger air 
mass. The overall refraction comes to about 0.5°, and calculating the path that light takes is hard, 
but you can find a bound on the answer by assuming a path that follows the surface of the Earth 
through this angle and then takes off on a straight line. What is the air mass ratio in this case? The 
real answer is somewhere between the two calculations. (The really real answer is a little bigger than 
either because the atmosphere is not isothermal and so the approximation p = pQe~^)^ is not exact.) 
Ans: ^ ^jRi^jlh = 32, +Re/h-^ 37. 
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13.35 Work in a thermodynamic system is calculated from dW = PdV. Assume an 
ideal gas, so that PV = nRT. (a) What is the total work, ^ dW, done around this cycle 
as the pressure increases at constant volume, then decreases at constant temperature, 
finally the volume decreases at constant pressure. 

(b) In the special case for which the changes in volume and pressure are very small, esti- 
mate from the graph approximately what to expect for the answer. Now do an expansion 
of the result of part (a) to see if it agrees with what you expect. Ans: ^ AP AV/2 

13.36 Verify the divergence theorem for the vector field 

F = axyzx + f3x'^z{l + y)y + -yxyz"^ z 
and for the volume {0 < x < a), {0 < y < b), {0 < z < c). 

13.37 Evaluate / F ■ dA over the curved surface of the hemisphere x'^ + y'^ + z"^ = and z > 0. 
The vector field is given by F = V x (^ayx + Pxy + '~fxyz) . Ans: {/3 — a)7ri?^ 

13.38 A vector field is given in cylindrical coordinates to be F = far'^z sin^ (p + cjy/Srz + Z'jzr cos^ 0. 
Verify the divergence theorem for this field for the region (0 < r < R), (0 < </> < 27r), {0 < z < h). 

13.39 For the function F{r,d) = r'^{A + BcosB + Ccos^O), compute the gradient and then the 
divergence of this gradient. For what values of the constants A, B, C, and (positive, negative, or zero) 
integer n is this result, V-VF, zero? These coordinates are spherical, and this combination divgrad 
is the Laplacian. 

Ans: In part, n = 2, C = -SA, B = 0. 

13.40 Repeat the preceding problem, but now interpret the coordinates as cylindrical (change 9 to (f)). 
And don't necessarily leave your answers in the first form that you find them. 

13.41 Evaluate the integral J F • dA over the surface of the hemisphere + + = 1 with z > 0. 

The vector field \s F = A{1 + x + y)x + B{1 + y + z)y + C{\ + z + x)z. You may choose to do this 
problem the hard way or the easy way. Or both. 
Ans: 7r(2v4 + 2fi + 5C)/3 

13.42 An electric field is known in cylindrical coordinates to he E = f{r)f, and the electric charge 
density is a function of r alone, p{r). They satisfy the Maxwell equation V ■ E = p/ Eq. If the charge 
density is given as p(r) = p^e'^/'^'^. Compute E. Demonstrate the behavior of E is for large r and 
for small r. 

13.43 Repeat the preceding problem, but now r is a spherical coordinate. 

13.44 Find a vector field F such that V ■ F = ax + /3?/ + 7 and V x F = 2. Next, find an infinite 
number of such fields. 

13.45 Gauss's law says that the total charge contained inside a surface \s eo ^ E ■ dA. For the electric 
field of problem 10.37, evaluate this integral over a sphere of radius Ri > R and centered at the origin. 

13.46 (a) In cylindrical coordinates, for what n does the vector field v = r"^ (p have curl equal to zero? 
Draw it. 

(b) Also, for the same closed path as in problem 13.17 and for all n, compute §v-dr. 
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13.47 Prove the identity Eq. (13.43). Write it out in index notation first. 

13.48 There an analog of Stokes' theorem for § M x B. This sort of integral comes up in computing 

the total force on the current in a circuit. Try multiplying (dot) the integral by a constant vector C . 
Then manipulate the result by standard methods and hope that in the end you have the same constant 
C • something. 

Ans: = J [(Vi?) ■ dA — (V ■ B) ■ dA~\ and the second term vanishes for magnetic fields. 

13.49 In the example (13.16) using Gauss's theorem, the term in 7 contributed zero to the surface 
integral (13.17). In the subsequent volume integral the same term vanished because of the properties 
of sm0cos(/>. But this term will vanish in the surface integral no matter what the function of (f) is in 
the (j) component of the vector F. How then is it always guaranteed to vanish in the volume integral? 

13.50 Interpret the vector field F from problem 13.37 as an electric field E, then use Gauss's law that 
^'enclosed = § E ■ dA to evaluate the charge enclosed within a sphere or radius R centered at the 
origin. 

13.51 Derive the identity Eq. (13.32) starting from the definition of a derivative and doing the same 
sort of manipulation that you use in deriving the ordinary product rule for differentiation. 

13.52 A right tetrahedron has three right triangular sides that meet in one vertex. 
Think of a corner chopped off of a cube. The sum of the squares of the areas of the 
three right triangles equals the square of the area of the fourth face. The results of 
problem 13.6 will be useful. 




